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' .'Abstract 

J-^ The length of minimal and maximal blocks equally distant on log-log scale versus fluctuation function considerably influences bias 
^ and variance of DFA. Through a number of extensive Monte Carlo simulations and different fractional Brownian motion/fractional 
Gaussian noise generators, we found the pair of minimal and maximal blocks that minimizes the sum of mean-squared error of 
estimated Hurst exponents for the series of length = 2*^,^ = 7, . . . , 15. Sensitivity of DFA to sort-range correlations was examined 
using ARFIMA(p, d, q) generator. Due to the bias of the estimator for anti-persistent processes, we narrowed down the range of 
(y;) Hurst exponent tol/2<_ff<l. 
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o 
o 

Introduction 

^ As of September, 2004, the two original paper s [T|2] on Detrended Fluctuation Analysis (DFA) had been cited by 470 
research articles, and the number is still growing]^ But still there is a need for improving the methodology, e.g., testing 
the limitations of DFA for various types of nonstationarities investigating its performance for signals with different 
type of correlations, with random spikes and missing segments, comparing it with other methods |4|5j , testing the effects 
of periodic (and quasi-periodic) trends in the estimation process |6l7j , studying the limitations of DFA for anti-persistent 
signals and the strategies to overcome them [8], its application to the wide class of multifractal series [9] or latest works 
on an impact of coarse-graining [10]. The method is also known as a Scaled Windowed Variance - Linear Detrended [TT] . 
Roughness Around the Root Mean Square Line [T^] and Residuals of Regression [13] . 

Kjj) Thorough examination of DFA as a type of linear regression detrended Scaled Windowed Variance method was carried 
• ^ out by [2] and later by [T3] Authors indicated sensitivity of DFA to exclusion of blocks of different size and after 
simulation stated that 

^ ,,. . . excluding windows of large sizes reduces variance but results in significant bias, while excluding windows of small 
sizes reduces bias in estimates nearly to zero for all signal lengths and all values of true H but the variance increases 
dramatically." 

In this study we will try to find the best block cuts by conducting exhaustive experiments with 10,000 replications using 
different types of fractional Brownian motion or fractional Gaussian noise generators. Final results will be verified under 
the presence of short-range dependence using ARIFMA(p, d, q) process. 
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2. Some preliminaries 



Self-similar real-valued stochastic process X — {X{t)}t£m. with Hurst exponent H > (H-ss) is defined as satisfying 

{X{at)}teM = {a^ X{t)}teM, for a > 0. Hyperbolically-decaying autocorrelation function j{k) of a stationary stochastic 
process {Xj}^q is nonsunimable (i.e. t(^) ~ defines asymptotically self-similar process if 

7(fc) cx fc2^-2^(fc), k^oo, 1/2<H<1, (1) 

where L{k) is a slowly-varying function, i.e. limj^oo L{tk) / L{k) = 1, and defines exactly self-similar process if 

7(/c) = l/2[(fc + - 2^" + {k- (2) 

If 7(A:) is diverging, one says that {Xt}ti£Z exhibits long-range dependence (LRD), strong dependence, has long memory 
or is a 1// noise '[T6\. LRD corresponds to the blow-up of the spectral density S{f) at the origin 

S{f)^cf-^", f^O, 0<H<1, (3) 

where 

^ oo 

^(/) = ^ E e"*^S(fc). (4) 

k — -'Oo 

Partitioning self-similar process {^iji^z into non-overlapping blocks of m sequential elements and computing average 
of these m elements 

i={t-l)m+l 

does not change the autocorrelation function |17j (contrary to "typical" stochastic processes — m increases and auto- 
correlation of {xj:"^^} decreases). 

The most widely-studied self-similar processes are fractional Gaussian noise (fGn) [TH] and autoregressive fractional 
integrated moving average processes (ARFIMA) jl5ll20j . In this study we used ARFIMA{0, d, 0), to generate H-ss. The 
general ARFIMA{p, d, q) process is defined as 

<f{L){l-LYXt^Q{L)et, (6) 

where L is the lag operator, Ct is white noise process and (i G K is the fractional differencing parameter such that 
|d| < 1/2. The process is covariance stationary if —1/2 < d < 1/2 and invertible for d > —1/2. For p,q — we have 
(1 — L^Xt — €t, and its Wold representation is given by 

oo 

X, = ^7^,L^ (7) 

3=0 

where ttq — I and 

Covariance function 7(fc) = EXoXk of ARFIMA{0, d, 0), -1/2 < d < 1/2 is given by [16] 

_ (-i)T(i-2rf) _ r(fc + d)r(i-2rf) 
^^''^"^ r(fc-d+i)r(i-fc-d) r(/fc-d+i)r(d)r(i-d) """'^l as/c^oo. (yj 

Hence, v4i?F/Af A(0, d, 0) has long memory if and only if < rf < 1/2. If we compare Eq. ([T|) and ^ we see that 
i? = d + 1/2. If {Xtjtez is a Gaussian ARFIMA{0, d, 0), < d < 1/2, then as n ^ oo n.-" ^ ^^(s), where 

{i?i/(s)}sgR is a fractional Brownian motion. 



3. Detrended Fluctuation Analysis 



Let Xi be the stationary series of compact support, where support is defined as a set of indices i with nonzero values 
Xi. The series is compact if Xi — for its small fraction only, interpreted then as having no value at this i. 

First the series is divided into [N/m\ non-overlapping^ logarithmically spaced blocks (windows, segments) of size 
m. Since N is often not a multiple of time scale to, a short part at the end of the series may remain. In order not to 

^ Overlapping blocks introduce correlations between estimates and should be abandoned. 
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disregard this part of the series, the same procedure is repeated starting from the opposite end, so 2N/m segments are 
obtained altogether. In this study we set N (and m) to a power of 2, so that N /m is directly integer and a final number 
of blocks. 

The advantage of DFA is that it is applied directly to nonstationary series. Hence, as we assumed Xi is stationary, the 
series must be integrated before analysis, calculating within each block partial sums Y{t) = -^i- During integration 

the sample mean X — N^^ T^^=i -^i '^^^ t)e subtracted (adjusted partial sums) Y{t) = Y^l^iiXi—X) but not compulsory 
because it will be eliminated by the later detrending. Within each fc = 1, 2, . . . , N/m block a least square line, + bkt, 
is fitted to the partial sums, and the sample variance of residuals is computed 

1 ™ / \^ 
F^{k,m) = ^ y((fc-l)m + t)-afe-5fctj , k = 1,2, N/m. (10) 

One can switch the fitting trend to quadratic, cubic, or higher order polynomials (called DFA2, DFA3, . . . , DFAr) [T] 
and by comparing the results for different r estimate the type of the polynomial trend in the time series [H] . Averaging 
Eq. ( 10 1 over all blocks gives the gth order fluctuation function 

1/9 

Fg{m) ~ - 



r , N/,n 



which is by construction defined only for m > r + 2. Eq. (Ill refers to multifractal detrended fluctuation analysis 
(MF-DFA) [22], thoroughly analyzed by Ref. [9 . In this study we focus only on the standard DFA, i.e. q — 2 

1/2 

I 1 I 

F(m) = ~ 



N/m ^ 

-Y^^[F\k,m))\ (12) 



N/ 

The variance of residuals is proportional to m^^ , where H is Hurst exponent [13], [T3]. Hence the fluctuation function 
described in Eq. (12 1 is proportional to 

F{m) (X m". (13) 

The scaling behavior of the fluctuation function is analyzed on log-log plots F{m) versus m - the slope of regression 
log(F(m)) = c + log TO is the Hurst exponent. As will be shown further DFA can reliably determine persistent signals 
(1/2 < i? < 1), and it becomes inaccurate for strongly anti-persistent processes {0 < H < 1/2) when H is close to zero. 
In such cases, a modified DFA technique has to be used. The easiest way is to apply — instead of single summation — 
double summation 

y(j)^^(r(i)-r). (i4) 



It leads to so called generalized fluctuation functions F{m) described by a scaling law with larger than in Eq. (13 1 
exponents 

F{m) (xm" ^ m"+^. (15) 

Comparing Eq. (13 1 and (15) we see that F{m)/m = F{m). Nevertheless double summation leads to quadratic trends 
in F{m). Hence if the average values were not removed in Eq. (14 1, at least the second order DFA should be applied 
to eliminate these artificial trends. Due to these inconveniences (and greater occurrence of persistent signals) we narrow 
down our analysis to persistent processes. 

The sum of the numbers X^ within each block k of size to is known as the box probability Pm{k) in the standard 
multifractal formalism for normalized series 

km 

p^{k)= Yl k = 1,2,..., N/m, (16) 

i={k-l)m+l 

and defines partition function with scaling exponent T(q) — qh{q) — 1 (here for q — 2) 

N/m 

Zg{m) = \Pmik)\'' oc TO«'^(«)-i = m^^-i. (17) 

k = l 

We can relate the scaling exponent from Eq. ( |17[ ) to Holder exponent a and singularity spectrum /(a) 

f{a) ^2{a- H) + l. (18) 
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4. Monte Carlo simulation 



In order to carry out the experiment we have chosen the following fractional Brownian motion generators: 

- Davies and Harte (known also as Wood and Chan, circulant matrix embedding method, exact) [25, [TB] , 

- Hosking method (recursive, exact, known also as Levinson method for Toeplitz matrices) j25], [26], 

- Choleski decomposition of the covariance matrix (exact) [26], [16], 
and fractional Gaussian noise generators: 

- Paxson method (approximate) |17j . 

- Beran method [?7] . 

- Durbin-Levinson (using Yule- Walker- type equations f^l]), [O], [55] , 

- ARFIMA(0, d, 0) (in the frequency domain using fast Fourier transform based on S-PLUS code written originally by 

my 

All fractional Gaussian noise series were cumulated before H estimation. Sensitivity of final results on short-range 
dependence was examined for ARFIMA(p, d, g) series: (0,c?, 1) with 9 = 0.5, {l,d,0) with = 0.5, (l,d, 1) with (j) = 
0.3, = 0.7, (1, d, 1) with (j) = -0.3, 6 = -0.7 and (1, d, 1) with = 0.7, d = 0.3. 

Because H is not known a priori, suggestions for block adjustment is based only on the length of the series. The 
optimization criterion has to be chosen to minimize both bias and variance for known N. Hence, a sum of mean-squared 
error (MSE) was chosen as the criterion of method reliability. Minimizing MSE accounts for square of bias and variance 
minimization [30' 

MSF.{H) = E{H - Hf = {EH - Hf +E{H -EHf = bia3^^ + variancefj (19) 

Proof 

MSE(^) = E{H - Hf = E{H^ + - 2HH) 
= E(i?2) + i/2 _ 2HEH 

= E(i?2) + ii-2 _ 2HEH + {EHf + {EHf - 2{EHf 
= {EHY + - 2HEH + E{H'^) + {EHf - 2{EHf 
= {EHf + - 2HEH + E{H^) + E{EHf - 2EHEH 
= {EHf + - 2HEH + E{H^ + {EHf - 2HEH) 
= {EH - Hf +E{H -EHf □ 

Let us consider MSE(iJ) as a function of a pair of minimal and maximal blocks {m~,m'^) chosen in the estimation 
process (Figure [ij . Assuming, that DFA is constructed on at least c* different equally distant on log- log scale blocks, we 
can define the set of all possible combinations of {m^,m'^) as 

C = {{m-,m+) = (2', 2"): u-l+l>c*Al = Zi, . . . , logj iV A u = ui, . . . , log2 iV}. (20) 




Fig. 1. An example of the fluctuation function F{m) and minimal and maximal blocks m , m+ (left and right dotted vertical lines). The 
slope of the solid line represents H, the slope of the dashed line - nominal H . 



4 



Table 1 

Number of elements of C for series of length N = 2"^ , . . . , 2^^ for different order of fitting polynomial trend r. For linear (r = 1) and quadratic 
(r = 2) polynomial trend li = 2. For cubic (r = 3) and 4th order polynomial fitting trend (r = A) li = 3. c* = 4. 



p 


\T OP 

i V — 










fitting trend 








linear or quadratic cubic 


or 4th order 


7 


128 


6 


3 


8 


256 


10 


6 


9 


512 


15 


10 


10 


1 024 


21 


15 


11 


2 048 


28 


21 


12 


4 096 


36 


28 


13 


8 192 


45 


36 


14 


16 384 


55 


45 


15 


32 768 


66 


55 



The number of all elements of C (presented in Table [T|) is then given by 

#C = a(a+l)/2, where a = log2iV + 2-?i -c*. (21) 

If the fitting trend is linear or quadratic (see Eq. ( [To| ) we set the shortest block to 2'^ = 4. For cubic and 4th order 
polynomial trend we set it to 2'^ = 8. The minimal number of different blocks is set to c* = 4. Let us also introduce the 
following function 

i?(m-,m+) = MSE(F)(m-,m+), (22) 



E 

Hen 



which idea is to describe the behavior of MSE(i7) for an a priori unknown Hurst parameter. Due to the strong bias of 
DFA for anti-persistent processes [5] we have considered in our simulation persistent processes only. Of course MSE(jy) 
is sensitive to length of series N and H, but to make the notation in Eq. (22 1 clear we omitted them. Hence, all values of 
d that will appear in our study are computed on H = {0.5, 0.6, 0.7, 0.8, 0.9}. Our goal is to find — through the number 
of computer simulations — the pair of (m^, m^), which will minimize d for long memory time series of length N 

(m~, m^) = arg min ?9(to~, rn'*'). (23) 

Before this step we tried approximate the number of the series iterations (replications), under which (m~,m+) would 
not change (within a small error). Hence, for a fixed number of replications of the fBm or fGn series of length N, we 
calculated (m~, m"*") and then looped this step 50 times, each time recording (m~, m"'")ea;, where ex indicates subsequent 
experiments. Then from these 50 optimal pairs of minimal and maximal blocks we calculated the mode and the frequency 
at which it occurred (in these 50 experiments) 

P{fh^ ,fn^) — max P ((m^, m+)ea;) , ea: = 1, . . . , 50, (24) 

(m^ ,m+)tEC 

and next we increased the number of replications. In Table [2] there is an example of the way P(to~, to"*") is determined 
for the series of length N = 128 iterated certain times. 



Table 2 

Example of the construction of P(m~,m"'") for TV = 128. Before the simulation is carried out we list C, which has six possible block 
combinations (m~,m+). After each of 50 experiments, (m~,m+) are sorted by •& in ascending order and P(m~,m"'") for the best (bold font) 
combination (m~,rh+) is computed. 



C after experiment 







1st 


2nd 


50th 


mode Pijh ,m+) 


(m , 


m+) 


(m~, m+) 


(m^, m+) 


(m~, m+) 


(m~, rh+) 


(4, 


128) 


(4, 32) 


(4, 64) 


(4, 32) 


(4, 32) 0.48 


(4, 


64) 


(4, 64) 


(4, 32) 


(8, 128) 




(4, 


32) 


(4, 128) 


(8, 64) 


(4, 128) 




(8, 


128) 


(8, 128) 


(4, 128) 


(4, 64) 




(8, 


64) 


(16, 128) 


(8, 128) 


(16, 128) 




(16, 


128) 


(8, 64) 


(16, 128) 


(8, 64) 
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Fig. 2. P(m~,m+) for A'' = 128 and 50 experiments against increasing iterations (by 10). While the first point is based on 250 fBm 
simulated series (10 iterations X 50 experiments X 5 Hurst exponents), the last one — on 2,500,000. Inner panel shows stability of mode (3 is 
(m~,m"'") = (4,32)). fBm generator: Davies and Harte. 

Figure [2] depicts the behavior of P{ih~ ,rh^) against the number of replications from 10 to 10,000 (increased by 10) 
and looped 50 times. Note that the last point on Figure |2] represents 2,500,000 generations oi fBm of length N = 128 
(10,000 iterations x 50 experiments x 5 Hurst parameters). We see that for 100 iterations P(to~,to+) is about 0.60, 
which means that only 60% of 50 carried experiments gave the same (to~, to"*"). For 1,000 iterations P(m~, to"*") reaches 
about 0.75 and for 10,000 is above 0.95. At this number of iterations there exists a small (less than 5%) error/risk that 
obtained (m^, m+) will change if the experiment is repeated. Hence, we will carry out our computations of i?(m^, m+) 
simulating 10,000 replications of LRD series for each oi H G {0.5, 0.6, 0.7, 0.8, 0.9}. For the purpose of a graphical 
presentation of the behavior of DFA H — {0.1,0.2, . . . ,0.9} has been chosen. Additionaly we will list up to three the 
best block combinations for different N and H. 

5. Results and discussion 

Table [3] presents the best three minimal and maximal block combinations (m~,m+) (notation (m^,m+) is reserved 
for #1) for 10,000 fBm/fGn independent paths generated with the use of different generators listed in sectionjl] Despite 
the differences in used generators, in five out of six the pair (4,32) minimizes i^(m~, m"*"), the second and third best 
pairs are (4,64) and (4,128) respectively. The only exception is the result obtained using Hosking generator, but the 
pairs (4, 32), (4, 64), (4, 128) are still preferable. We also see that the (m^, m+) do not change with the increase of the 
number of observations N. 

'd{m~ ,m~^) for different generators for the same N and rank {^1,^2,^3) are very close to each another, which 
indicates that DFA behaved in similar way on the simulated processes. Hence, to picture bias, standard deviation and 
root mean-squared error of H we have chosen Davies and Harte generator of fractional Brownian motion. The generator 
was one of the fastest and that is why we simulated up to 32,768 observations. Figure |3] depicts boxplots of the bias 
of the estimated Hurst parameter H = 0.1, . . . ,0.9 on 10,000 fbm series of length N = 32768 for all possible block 
combinations. In Table [4] we listed bias, standard deviation and RMSE for the best three (TO~,m+). To make it more 
readable we presented the output on Figures [5](bias),|6](std. dev.) and [t] (RMSE) as a function of H, log2 N. On lower 
panels of the Figures |5]|7] we plot cubic splines with their contours for H x logj N. For the purpose of a visual presentation 
of the behavior of DFA we extended the set of available Hurst parameters to H — {0.1, . . . ,0.9} bearing in mind that 
d{rh~ ,rh~^) was computed on H = {0.5, . . . , 0.9}. 

Let us give a short description of the behavior of the estimator. Let's consider (m^, m+) — (4, 32). For the series of 
length N = 2^°, . . . , 2^^ with nominal H = 1/2 bias is about 0.004 - 0.005, while for such a span of N standard deviation 
decreases from 0.029 to 0.005 and RMSE decreases from 0.029 to 0.007 (Table [i]). With the increase of H from 0.5 to 0.9 
bias changes from 0.004 to —0.014 for N = 1024 and from 0.005 to —0.010 for N — 32768. Standard deviation increases 
from 0.005 for 7? = 1/2 to 0.007 for H = 0.9 for longest available series. The behavior of DFA for H = 0.1, 0.2, 0.3, 0.4, 
which is not presented in tables, can be seen on Figures [5] - [7] 

Let us also compare obtained results using the optimal pair of minimal and maximal blocks: (m^,m+) = (4,32) and 
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the pairs suggested by Ref. HM- (8, 256) for N = 1024; (16, 256) for N ^ 8192; (64, 1024) for N = 32768. In our simulation 
these pairs are respectively: #8 out of 21, #13 out of 45 and ^^32 out of 66 possible combinations. We have compared 
bias and standard deviation for such scales in Table |5] for N — 1024 and N — 8192 and in Table |6] for N — 32768. Let 
us briefly describe the results. For the series of length N — 1024 bias for (m^, m+) = (4, 32) is lower from 1.7 up to 4.7 
times than in the (8, 256) combination and standard deviation is lower about 1.5 times. With the increase of the length 
of simulated paths {N — 8192) the difference in bias is getting smaller — 4.7 times bigger for H = 0.6 and equal for 
H = 0.9, but the difference in standard deviation increases to two times. For the longest series available {N — 32768) 
bias is lower (about a half) in favor of the (64, 1024) combination (except for H — 0.6) but standard deviation rose 
almost up to 4 times than in the optimal pair. That is why i9(64, 1024) was 32nd out of 66 possible combinations. 



Table 3 

The best three block combinations {•m~,m~^) and ^{m~ ,m~^) for fBm/fGn series of length TV = 2'^, . . . , and H = 0.5, 0.6, . . . , 0.9. With 
X we marked not simulated — due to complexity of the algorithm — series of length A'^. Generators: hos - Hosking, dh - Davies and Harte, 
chol - Cholcsky decomposition, pax — Paxson, ber - Beran, dl - Durbin-Levinson. 



N 














(m 


, m+) 
















m , rh'^ 


) 






hos 




dh 


chol 


pax 


ber 




dl 




dh 








dl 


128 


#1 


(4, 


64) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


0.0465 


0.0493 


0.0495 


0.0484 


0.0481 


0.0483 




#2 




32) 


(4, 


64) 




64^ 


(d 




(d 


64") 


(4, 


64) 


0.0479 


0.0499 


0.0501 


0.0493 


0.0482 


0.0486 




#3 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


0.0515 


0.0567 


0.0560 


0.0573 


0.0545 


0.0558 


256 


#1 


(4, 


64) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


0.0236 


0.0252 


0.0248 


0.0246 


0.0245 


0.0247 




#2 




32) 


(4, 


64) 




64) 


(d 


64) 


(d 


64) 


(4, 


64) 


0.0243 


0.0254 


0.0251 


0.0247 


0.0250 


0.0248 




#3 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


0.0261 


0.0289 


0.0285 


0.0287 


0.0286 


0.0284 


512 


#1 


(4, 


64) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


0.0122 


0.0127 


0.0125 


0.0124 


0.0125 


0.0126 




#2 


(4, 


32) 


(4, 


64) 


(4, 


64) 


(4, 


64) 


(4, 


64) 


(4, 


64) 


0.0128 


0.0130 


0.0127 


0.0127 


0.0127 


0.0129 




#3 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


0.0135 


0.0146 


0.0145 


0.0146 


0.0145 


0.0146 


1024 


#1 


(4, 


64) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


0.0066 


0.0065 


0.0065 


0.0064 


0.0064 


0.0064 




#2 


(4, 


32) 


(4, 


64) 


(4, 


64) 


(4, 


64) 


(4, 


64) 


(4, 


64) 


0.0068 


0.0067 


0.0067 


0.0065 


0.0066 


0.0065 




#3 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


0.0071 


0.0077 


0.0075 


0.0075 


0.0075 


0.0074 


2048 


#1 


(4, 


64) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


0.0036 


0.0033 


0.0033 


0.0033 


0.0033 


0.0033 




#2 


(4, 


32) 


(4, 


64) 


(4, 


64) 


(4, 


64) 


(4, 


64) 


(4, 


64) 


0.0037 


0.0035 


0.0034 


0.0034 


0.0034 


0.0035 




#3 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


128) 


0.0039 


0.0040 


0.0039 


0.0039 


0.0039 


0.0039 


4096 


#1 


(4, 


64) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 


32) 


(4, 
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H=0.1 H=0,2 M=0.3 




Fig. 3. Boxplots for 10,000 fBm paths of length A'^ = 32768 generated using Davies and Harte exact method. On X axis — different block 
combinations, starting from 4 obs. in the shortest block and N in the longest one. Longest blocks are cut first till at least four blocks are left. 
On Y axis — deviation of the nominal value. 



K = 12a, opt mQl blocks (-1,32: N = 255. optimal blocks (4,32) N=512, optimal blocks (4,32) 




23456759 1 2 3 4 5 6 7 8 9 1234567E9 



N=8192. optimal blocks (4,32) N=16384, optimal blocks (4,32) N=32768, optimal blocks (4,32) 




10.H 10.H 10»H 



Fig. 4. The optimal pair of minimal and maximal blocks (m ,m+) = (4,32) and boxplots of the bias of the estimated Hurst parameter 
H = 0.1,..., 0.9 for 10,000 fBm series of length N = 2"^ , . . . ,2^^ simulated by Davies and Harte method. Note: for if = 0.1, ... , 0.4 the pair 
(4,32) may not be optimal. 
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08 




Fig. 5. Upper left: scatter plot of bias against H with fitted lines for different length of the series N = 2^, . . . , 2^^. Upper right; scatter plot 
of bias against log2 N with fitted lines for different Hurst parameters, H = O.l, . . . , 0.9. Lower right: cubic spline fit of bias with H X log2 N. 
Lower left: contour of cubic spline fitted data. 10,000 fBm paths generated with Davies and Harte algorithm, rh~ = 4, m+ = 32. Note that 
for Ti. = {0.1, 0.2, 0.3, 0.4} (4, 32) may not be the optimal pair — the values are shown just for a visual presentation. 




Fig. 6. Same as Figure [5] but instead of bias — standard deviation. 
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Fig. 7. Same as Figure [s] but instead of bias — root mean-squared error. 

Table 4 

Bias, standard deviation and root moan square error for the best three block combinations. 10,000 estimates of H for different lengths of N , 
minimal m~ and maximal m+ blocks. Process: /Srn, generator: Davies and Harte. 
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N = 8<92. blocks (16,256) 



N = S192. ootimol blocks (4,32) 




Fig, 8, Comparison of boxplots of the bias of the estimated Hurst parameter for the the optimal pair of minimal and maximal blocks (4, 32) 
and for (64, 1024), 10,000 paths of length N = 8192 of fBm simulated by Davies and Harte method. 

Table 5 

Comparison of bias and standard deviation for the (m~ , rh+) = (4, 32) and (8, 256), (16, 256) block combinations for the series of N = 1024 
and N = 8192 respectively. 



H N = 1024 N = 8192 





Bias 






Std. dev. 






Bias 






Std. dev. 






(8,256) 
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Note: (a), (6), (a/b) are rounded numbers. 



Table 6 

Comparison of bias and standard deviation for the (4,32) and (64, 1024) block combinations for the series of length N = 32768, 

H N = 32768 



Bias Std. dev. 





(64, 1024) 
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(64,1024) (4,32) 




(a) 


(6) {a/b) 


(a) (b) {a/b) 
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0.024 0.006 3.8 


0.8 
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-0.004 


-0,010 0.4 


0.027 0.007 3.9 



5,1, The polynomial trend fit 

Tables and figures for this and the following subsection were moved to appendix \K\ 

We have also tried to investigate the impact of different orders of polynomial fit on the properties of DFA, Hence, 



we have chosen hnear, quadratic, cubic and fourth order polynomial fit. The results are presented in Table |A,1| Our 
optimization criterion d varies slightly for long series. One can use 4th order polynomial fit for the series of length of 
N = 32768 and (m", m+) = (16, 256) where bias is 0.006 for H = 1/2 and 0.000 for H = 0.8, 0.9 with standard deviation 
only approximately 0.002 grater than in the case of using the linear fit, but should not use r = 4 for e.g. N=1024, where 
bias for H = 1/2 is 0.026 comparing to 0.004 for r = 1 — six time higher. For the simplicity of empirical analysis we 
suggest to use the simplest case which is linear trend fit. Note that {m~ ,7ti~^) is the same for different iV only for that 
kind of fit, what is more convenient to carry out empirical analysis. 
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5.2. Impact of short memory 



— ARFIM A (p,d,q) 



We have examined DFA on the simulated ARIFMA(p, d, q) series in the separate subsection due to the following 
reasons 

- The asymptotic expansion of ARIFMA(0,d,0) processes, which may influence final results [23]. 

- Ability to verify the impact of additional short-range correlations, i.e., p,q > 0, on the blocks adjustment. 

We have simulated ARFIMA(p, d, q) in the frequency domain using fast Fourier transform based on S-PLUS code 
written originally by [27j and then checked sensitivity of final results on short-range dependence for ARFIMA(p, d, q) 
series: (0,d,l) with 9 = 0.5, (l,rf,0) with <j) = 0.5, 1) with 0.3,6* = 0.7, 1) with = -0.3,61 = -0.7 and 



1) with (j) = 0.7,6* = 0.3. Summary results for the best three block combinations are in Table A. 2 
Let us describe obtained results. DFA behaves differently on autoregressive fractional integrated moving average 
process than on fractional Gaussian noise or fractional Brownian motion. For ARFIMA(0, d, 0) with up to = 512 
observations, among the best three block combinations there is (4,32) (the pair (4,64) minimizes -d). In Table A. 3 we 
see that — contrary to bias — standard deviation for different Hurst parameters is almost identical. Except for H — 0.5, 
bias is much larger for ARFIMA(0, d, 0) than for Davies and Harte generator and that is why (TO^,m+) — (4,32) is no 
longer valid for such process. 

In Table [Ar4] we have presented the sensitivity of DFA on the different values of autoregressive (j) and moving average 9 



parameters for the same block combination (4,32). Introducing MA part leads to strong negative bias (see Figure A. 2 1 
while AR causes strong positive bias (Figure 



than (4, 32) optimal pairs listed in Table A. 5 



A. 3 1 (variance of the estimator remains stable). This results in the different 
Although we managed to decrease the bias, standard deviation rose, which 
makes DFA incapable to estimate Hurst parameter precisely in terms of the presence of short memory. 

DFA "prefers" much larger scales (m~, m"*") under the presence of short memory, but for the same length of the series 
these scales are different (Table [A3| . After the joint introduction of AR and MA parts we still observe strong bias and 



steady variance. For = 0.7 > 9 = 0.3 bias is positive and if (j) ~ 0.3 < 9 ^ 0.7 — negative (Figures A. 4 A.6 1. Negative 



(f) = —0.3 and 9 = —0.7 parameters (Figure A. 5) lead to greater bias than in the case of lack of short memory, but the 
effect is not as strong as for positive values of (f) and 9 (and smaller scales are preferable). 

Due to the complexity of the algorithm, we have restricted that part of our analysis to six different {p, d, q) only (three 
for (1, d, 1)). To obtain the full picture of the behavior of DFA on ARFIMA(0, d, 1), (1, d, 0) and (1, d, 1) we should carry 
out our simulations for </>, 6* = 0, ±.1, ... , ±0.9, d = 0, ±0.1, . . . , ±0.4 and = 2^, . . . , 2^. It took about 24 hours to 
simulate 10,000 certain ARFIMA(1, d, 1) series of length N — 8192 with nine different values of d parameter. Hence, 
such extended and precise analysis (in terms of the number of replications) would be very time-consuming. Nevertheless 
we have tried to find best block combinations on available six simulated ARFIMA(p, d, g) processes, searching for the 
pair (to^, TO+) that minimizes the sum of MSE for H = 0.5, . . . , 0.9 and these six processes altogether, bearing in mind 
that the results will be preliminary. We have listed them in Table |A.6[ Although we found the combination for the series 
of length N , we see that the variance of DFA is too large (standard deviation for 8192 observations is about 0.04-0.06 
with bias varying form —0.09 to 0.029). 

6. Conclusions 

Blocks adjustment may significantly improve precision of DFA for persistent processes. These improvements are pre- 
sented in Tables |5] |6] and depicted on Figures [3] |4] and |8] For the sake of simplicity of analysis we recommend using 
linear trend fit (Table [ATT] ). Through extensive simulations we have shown that the optimal pair of minimal and maximal 
blocks is (4,32), which may reduce standard deviation even up to four times for the series of length N — 32768. Second 
and third results are pairs (4,64), (4, 128) respectively. The results are robust on the length of the series and the type 
of generators of fractional Brownian motion or fractional Gaussian noise. The exception is ARFIMA(0, d, 0) — blocks 



adjustment also improved the quality of the estimator, but indicated different block combinations as optimal (Table A.2 1. 

DFA is very sensitive to the presence of short-range correlations in the series. The bias and variance are reduced at 
the expense of the other — exclusion of large-sized blocks reduces variance but results in greater bias, contrary to cuts 
of small blocks — but the effect is so strong that it disables precise estimation of the long memory parameter. In such 
cases data must be filtered or DFA should be replaced e.g., with Global log-periodogram 16J. 
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Appendix A. Tables and figures for sections |5.1| and |5.2| 



Table A.l 

Impact of order r = 1,2, ...,4 of the polynomial trend fit on bias, standard deviation and root mean-squared error for the best block 
combination {m~,m~^). Generator: Davies and Ifarte method for 10,000 paths of fractional Brownian motion. 
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041 





043 





046 





033 





037 





041 





044 





047 




3 


(8, 


256) 





0072 





013 





009 





005 





002 





000 





031 





034 





038 





040 





043 





034 





036 





038 





040 





043 




4 


(8, 


512) 





0085 





026 





021 





017 





013 





Oil 





031 





034 





038 





039 





042 





041 





040 





041 





041 





044 


2048 


1 


(4, 


32) 





0033 





004 


-0 


002 


-0 


007 


"0 


Oil 


-0 


012 





021 





023 





025 





026 





028 





021 





023 





026 





028 





030 




2 


(8, 


128) 





0042 





002 


-0 


003 


-0 


005 


-0 


007 


-0 


008 





024 





026 





028 





031 





032 





024 





026 





029 





032 





033 




3 


(8, 


256) 





0040 





016 





010 





007 





005 





004 





022 





024 





027 





029 





030 





027 





027 





028 





029 





030 




4 


(16, 


256) 





0051 





005 





001 


-0 


001 


-0 


003 


-0 


003 





026 





029 





032 





035 





036 





027 





029 





032 





035 





036 


4096 


1 


(4, 


32) 





0018 





005 


-0 


003 


-0 


007 


-0 


009 


-0 


Oil 





015 





016 





017 





019 





020 





015 





016 





019 





021 





023 




2 


(8, 


128) 





0021 





002 


-0 


002 


-0 


004 


-0 


005 


-0 


007 





017 





018 





020 





022 





023 





017 





019 





021 





022 





024 




3 


(8, 


256) 





0023 





016 





Oil 





009 





007 





005 





016 





017 





019 





020 





021 





022 





021 





021 





022 





022 




4 


(16, 


256) 





0026 





005 





002 





001 





000 


-0 


002 





019 





021 





023 





024 





026 





019 





021 





023 





024 





026 


8192 


1 


(4, 


32) 





0010 





004 


-0 


002 


-0 


006 


-0 


009 


-0 


Oil 





010 





Oil 





012 





013 





014 





Oil 





Oil 





014 





016 





017 




2 


(8, 


128) 





0011 





002 


-0 


001 


-0 


004 


-0 


005 


-0 


006 





012 





013 





014 





015 





016 





012 





013 





015 





016 





017 




3 


(8, 


512) 





0013 





Oil 





008 





006 





004 





003 





012 





013 





014 





016 





017 





017 





016 





016 





016 





017 




4 


(16, 


256) 





0013 





006 





003 





001 





000 


-0 


001 





013 





015 





016 





017 





019 





014 





015 





016 





017 





019 


16384 


1 


(4, 


32) 





0006 





004 


-0 


002 


-0 


006 


-0 


009 


-0 


010 





007 





008 





009 





009 





010 





009 





008 





Oil 





013 





014 




2 


(8, 


128) 





0006 





002 


-0 


001 


-0 


003 


-0 


005 


-0 


005 





008 





009 





010 





Oil 





012 





009 





009 





Oil 





012 





013 




3 


(8, 


1024) 





0008 





008 





006 





004 





003 





002 





010 





Oil 





012 





013 





014 





013 





012 





013 





013 





014 




4 


(16, 


256) 





0007 





006 





003 





001 





000 





000 





009 





Oil 





Oil 





012 





013 





Oil 





Oil 





012 





012 





013 


32768 


1 


(4, 


32) 





0004 





005 


-0 


002 


-0 


006 


-0 


009 


-0 


010 





005 





006 





006 





007 





007 





007 





006 





009 





Oil 





012 




2 


(8, 


128) 





0003 





003 


-0 


001 


-0 


003 


-0 


004 


-0 


005 





006 





007 





007 





008 





008 





006 





007 





008 





009 





010 




3 


(16, 


256) 





0004 





000 


-0 


002 


-0 


003 


-0 


004 


-0 


004 





007 





008 





009 





010 





010 





007 





008 





010 





010 





Oil 




4 


(16, 


256) 





0004 





006 





003 





001 





000 





000 





007 





007 





008 





009 





009 





009 





008 





008 





009 





009 
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Table A.2 

The best three block combinations (m-,m+) for 10,000 ARFIMA(p, d, g) series of length N = 2^, . . . , and if = 0.5, . . . , 0.9. 

AT (m~,m+) 



(0,d,0) (0,d, 1) (l,<i,0) (l,d,l) 1) (1,^,1) 

61 = 0.5 <p = 0.5 <t> = 0.3 <p = -0..3 <p = 0.7 

e = 0.7 6» = -0.7 = 0.3 



128 


#1 


(4, 


64) 


(16, 


128) 


(8, 


12) 


(16, 


128) 


(4, 


128) 


(8, 


128) 




#2 


(4, 


32) 


(8, 


128) 


(4, 


128) 


(4, 


128) 


(4, 


64) 


(16, 


128) 




#3 


(4, 


128) 


(4, 


128) 


(16, 


128) 


(8, 


128) 


(4, 


32) 


(4, 


128) 


256 


#1 


(4, 


64) 


(16, 


256) 


(8, 


256) 


(32, 


256) 


(4, 


128) 


(16, 


256) 




#2 


(4, 


32) 


(32, 


256) 


(16, 


256) 


(16, 


256) 


(4, 


64) 


(8, 


256) 




#3 


(4, 


128) 


(16, 


128) 


(8, 


128) 


(8, 


256) 


(4, 


256) 


(32, 


256) 


512 


#1 


(4, 


64) 


(32, 


512) 


(8, 


512) 


(32, 


512) 


(4, 


128) 


(16, 


512) 




#2 


(4, 


128) 


(16, 


512) 


(16, 


512) 


(64, 


512) 


(4, 


256) 


(32, 


512) 




#3 


(4, 


32) 


(32, 


256) 


(16, 


256) 


(32, 


256) 


(4, 


64) 


(8, 


512) 


1024 


#1 


(4, 


128) 


(32, 


1024) 


(8, 


1024) 


(64, 


1024) 


(4, 


256) 


(16, 


1024) 




#2 


(4, 


256) 


(32, 


512) 


(16, 


512) 


(32, 


1024) 


(4, 


512) 


(32, 


1024) 




#3 


(4, 


64) 


(16, 


1024) 


(16, 


1024) 


(64, 


512) 


(4, 


128) 


(32, 


512) 


2048 


#1 


(4, 


256) 


(32, 


1024) 


(16, 


1024) 


(64, 


1024) 


(4, 


512) 


(16, 


2048) 




#2 


(4, 


128) 


(32, 


2048) 


(8, 


2048) 


(64, 


2048) 


(8, 


128) 


(32, 


1024) 




#3 


(4, 


512) 


(64, 


1024) 


(16, 


2048) 


(32, 


2048) 


(4, 


256) 


(32, 


2048) 


4096 


#1 


(4, 


512) 


(64, 


1024) 


(16, 


2048) 


(64, 


2048) 


(8, 


128) 


(32, 


2048) 




#2 


(4, 


256) 


(64, 


2048) 


(16, 


1024) 


(64, 


4096) 


(8, 


256) 


(16, 


4096) 




#3 


(8, 


256) 


(32, 


2048) 


(16, 


4096) 


(64, 


1024) 


(8, 


64) 


(32, 


4096) 


8192 


#1 


(8, 


512) 


(64, 


2048) 


(32, 


1024) 


(128 


2048) 


(8, 


256) 


(32, 


4096) 




#2 


(4, 


512) 


(64, 


1024) 


(16, 


4096) 


(64, 


4096) 


(8, 


128) 


(32, 


2048) 




#3 


(4, 


1024) 


(64, 


4096) 


(32, 


2048) 


(64, 


2048) 


(8 


512) 


(32, 


8192) 
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Table A.3 

Comparison of the behavior of DFA for the pair (4, 32) on fractional Brownian motion simulated using Davies and Harte exact method (dli) 
and on ARFIMA(0, d, 0). 



N Numiual 11 









0.5 




0.6 




0.7 




0.8 




0.9 


0.5 


0.6 




0.7 




0.8 




0.9 


Bias 


















Std. dev. 














128 




dh 


-0.005 


-0, 


.013 


-0 


.024 


-0, 


.028 


-0, 


.037 


0.082 


0.089 


0, 


,096 


0, 


,102 


0, 


,110 




ARFIMA(0, 


d, 0) 


-0.003 


-0, 


.030 


-0, 


.053 


-0, 


,072 


-0 


.083 


0.081 


0.089 


0, 


,095 


0, 


,103 


0, 


,110 


256 




dh 


0.000 


-0, 


.009 


-0, 


.015 


-0, 


.020 


-0, 


.025 


0.058 


0.064 


0, 


,069 


0, 


,074 


0, 


,079 




ARFIMA(0, 


d,0) 


0.000 


-0, 


.026 


-0 


.046 


-0, 


.061 


-0 


.074 


0.058 


0.062 


0, 


,067 


0, 


,073 


0, 


,077 


512 




dh 


0.002 


-0, 


.005 


-0, 


.011 


-0, 


,014 


-0, 


,018 


0.041 


0.045 


0, 


,049 


0, 


,053 


0, 


,056 




ARFIMA(0, 


d,0) 


0.002 


-0, 


.023 


-0, 


.043 


-0, 


,059 


-0, 


.069 


0.041 


0.045 


0, 


,049 


0, 


,051 


0, 


,055 


1024 




dh 


0.004 


-0, 


.004 


-0 


.009 


-0, 


,012 


-0 


.014 


0.029 


0.032 


0, 


,035 


0, 


,037 


0, 


,039 




ARFIMA(0, 


d,0) 


0.004 


-0, 


.021 


-0, 


.042 


-0, 


,056 


-0, 


,067 


0.029 


0.031 


0, 


,034 


0, 


,036 


0, 


,039 


2048 




dh 


0.004 


-0, 


.002 


-0, 


.007 


-0, 


,011 


-0, 


.012 


0.021 


0.023 


0, 


,025 


0, 


,026 


0, 


,028 




ARFIMA(0, 


d,0) 


0.004 


-0, 


.021 


-0 


.041 


-0, 


,056 


-0, 


.066 


0.021 


0.023 


0, 


,024 


0, 


,026 


0, 


,028 


4096 




dh 


0.005 


-0, 


.003 


-0, 


.007 


-0, 


,009 


-0, 


,011 


0.015 


0.016 


0, 


,017 


0, 


,019 


0, 


,020 




ARFIMA(0, 


d,0) 


0.004 


-0, 


.021 


-0, 


.040 


-0, 


,055 


-0, 


.065 


0.015 


0.016 


0, 


,017 


0, 


,018 


0, 


,019 


8192 




dh 


0.004 


-0, 


.002 


-0 


.006 


-0, 


,009 


-0 


.011 


0.010 


0.011 


0, 


,012 


0, 


,013 


0, 


,014 




ARFIMA(0, 


d, 0) 


0.004 


-0, 


.021 


-0, 


.040 


-0, 


,055 


-0, 


,065 


0.010 


0.011 


0, 


,012 


0, 


,013 


0, 


,014 


16384 




dh 


0.004 


-0, 


.002 


-0, 


.006 


-0, 


,009 


-0, 


.010 


0.007 


0.008 


0, 


,009 


0, 


,009 


0, 


,010 




ARFIMA(0, 


d,0) 


0.004 


-0, 


.021 


-0, 


.040 


-0, 


,055 


-0, 


.065 


0.007 


0.008 


0, 


,009 


0, 


,009 


0, 


,010 



Tabic A.l 

Example of the behavior of DFA for the jjair (1..32) on fDrn simulated by Davies and llartc' exact method (dh) and on ARFIMA(p, d, q) of 
length N = 8192. 



AT = 8192, (m" , m" 


^) = (4,32) 


Nominal H 






















0.5 


0.6 


0.7 


0.8 


0.9 


0.5 


0.6 


0.7 


0.8 


0.9 






Bias 










Std. dev. 








dh 




0.004 


-0.002 


-0.006 


-0.009 


-0.011 


0.010 


0.011 


0.012 


0.013 


0.014 


ARFIMA(0, d, 0) 




0,004 


-0.021 ■ 


-0.040 ■ 


-0.055 


-0.065 


0.010 


0.011 


0.012 


0.013 


0.014 


ARFIMA(0, d, 1) e 


= 0.5 


-0.276 


-0.323 


-0.361 ■ 


-0.389 


-0.409 


0.007 


0.008 


0.009 


0.010 


0.012 


ARFIMA(l,d,0) <f> 


i = 0.5 


0.311 


0.287 


0.266 


0.246 


0.228 


0.012 


0.013 


0.013 


0.014 


0.015 


ARFIMA(l,d, 1) rj) 


i = 0.3, 6» = 0.7 


-0.255 


-0.306 


-0.350 ■ 


-0.386 


-0.414 


0.006 


0.007 


0.008 


0.010 


0.011 


ARFIMA(l,(i, 1) 4> 


. = -0.3,61 = -0.7 


0.105 


0.078 


0.055 


0.036 


0.020 


0.011 


0.012 


0.012 


0.013 


0.014 


ARFIMA(l,(i, 1) 4> 


. = 0.7,6» = 0.3 


0.343 


0.328 


0.314 


0.301 


0.289 


0.013 


0.014 


0.014 


0.015 


0.015 



Table A.5 

Bias and standard deviation for DFA based on optimal blocks {m^ ,m,~) for fBm simulated by Davies and Harte exact method (dh) and 
ARFIMA(p, d, q) of length N = 8192. 



N = 8192 






(m , m ) 


Nominal H 


























0.5 


0.6 


0.7 


0.8 


0.9 


0.5 


0.6 


0.7 


0.8 


0.9 










Bias 










Std. dev. 








dh 






(4,32) 


0.004 


-0.002 


-0.006 


-0.009 


-0.011 


0.010 


0.011 


0.012 


0.013 


0.014 


ARFIMA(0, d, 0) 






(8,512) 


-0.009 


-0.017 


-0.021 


-0.023 


-0.024 


0.016 


0.018 


0.020 


0.021 


0.023 


ARFIMA(0, d, 1) e 


= 0.5 




(64, 2048) 


-0.058 


-0.045 ■ 


-0.036 


-0.030 


-0.027 


0.040 


0.045 


0.050 


0.054 


0.058 


ARFIMA(l,d,0) (j> 


= 0.5 




(32, 1024) 


0.032 


0.023 


0.015 


0.010 


0.006 


0.030 


0.033 


0.036 


0.039 


0.042 


ARFIMA(l,d, 1) (j) 


= 0.3,6* 


= 0.7 


(128, 2048) 


-0.068 


-0.052 


-0.043 


-0.033 


-0.031 


0.051 


0.059 


0.065 


0.072 


0.076 


ARFIMA(l,d, 1) </) 


= -0.3, 


e = -0.7 


(8,256) 


0.019 


0.007 


-0.001 


-0.007 


-0.011 


0.015 


0.017 


0.018 


0.020 


0.021 


ARFIMA(l,(i, 1) (t> 


= 0.7, 6» 


= 0.3 


(32, 4096) 


0.037 


0.026 


0.018 


0.011 


0.006 


0.037 


0.041 


0.046 


0.048 


0.051 
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Table A.6 

Bias and standard deviation for DFA based on optimal blocks (ni^ , fh^) for available ARFIMA(p, d, q) of length N = 2'^ , . . . , 2^^. 

{■m~,rh~) Nominal -f/ 

05 06 07 08 09 05 06 07 08 09 
Bias Std. dev. 



N = 128 



ARFIMA 





d, 


0) 






(8,128) 


-0 


037 


-0 


053 


-0 


068 


-0 


077 


-0 


087 





103 





115 





125 





136 





145 


ARFIMA 





d, 


1) 


e = 0.5 




(8,128) 


-0 


242 


-0 


263 


-0 


275 


-0 


276 


-0 


278 





080 





096 





111 





125 





138 


ARFIMA 


1 


d, 


0) 


<t> = 0.5 




(8,128) 





124 





101 





079 





059 





043 





114 





123 





131 





142 





150 


ARFIMA 


1 


d, 


1) 


4> = 0.3.6* 


= 0.7 


(8,128) 


-0 


275 


-0 


305 


-0 


326 


-0 


340 


-0 


345 





070 





087 





102 





115 





130 


ARFIMA 


1 


d. 


1) 


(f> = -0.3 


e = -0.7 


(8,128) 





000 


-0 


020 


-0 


039 


-0 


053 


-0 


061 





106 





117 





126 





136 





144 


ARFIMA 


1 


d, 


1) 


<t> = 0.7.9 


= 0.3 


(8, 128) 





199 





177 





156 





138 





114 





119 





129 





137 





145 





149 



N = 256 



ARFIMA 





d, 


0) 






(16,256) 


-0 


035 


-0 


043 


-0 


054 


-0 


062 


-0 


071 





101 





113 





124 





135 





146 


ARFIMA 





d, 


1) 


e = 0.5 




(16,256) 


-0 


185 


-0 


188 


-0 


185 


-0 


181 


-0 


170 





086 





103 





117 





130 





141 


ARFIMA 


1 


d, 


0) 


<t> = 0.5 




(16,256) 





060 





038 





021 





005 


-0 


009 





109 





119 





130 





139 





146 


ARFIMA 


1 


d, 


1) 


= 0.3. e 


= 0.7 


(16,256) 


-0 


234 


-0 


246 


-0 


248 


-0 


246 


-0 


238 





079 





096 





111 





125 





138 


ARFIMA 


1 


d, 


1) 


</> = -0.3 


e = -0.7 


(16,256) 


-0 


016 


-0 


032 


-0 


041 


-0 


051 


-0 


061 





103 





115 





123 





136 





146 


ARFIMA 


1 


d, 


1) 


</> = 0.7.61 


= 0.3 


(16,256) 





121 





099 





080 





061 





045 





114 





121 





131 





140 





149 



N = 512 



ARFIMA 





d, 


0) 






(16,512) 


-0 


028 


-0 


038 


-0 


044 


-0 


050 


-0 


056 





079 





087 





096 





104 





110 


ARFIMA 





d, 


1) 


6» = 05 




(16,512) 


-0 


158 


-0 


154 


-0 


150 


-0 


141 


-0 


131 





071 





083 





092 





102 





109 


ARFIMA 


1 


d, 


0) 


= 0.5 




(16,512) 





048 





030 





012 





001 


-0 


010 





081 





089 





098 





105 





114 


ARFIMA 


1 


d, 


1) 


4> = 0.3,6* 


= 0.7 


(16,512) 


-0 


204 


-0 


209 


-0 


206 


-0 


199 


-0 


191 





067 





079 





090 





099 





108 


ARFIMA 


1 


d, 


1) 


= -0.3, 


e = -0.7 


(16,512) 


-0 


013 


-0 


025 


-0 


036 


-0 


042 


-0 


051 





079 





087 





096 





104 





111 


ARFIMA 


1 


d, 


1) 


(t> = 0.7,6 


= 0.3 


(16,512) 





099 





082 





063 





048 





031 





084 
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Fig. A.l. Boxplots for 10,000 ARFIMA(0,d,0) series of length = 8192. On X axis — different block combinations, starting from 4 obs. in 
the shortest block and A^ in the longest one. Longest blocks are cut first till at least four blocks are left. On Y axis — deviation of the nominal 
value. 
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Fig. A.2. Same as Figure [aTT] but ARFIMA(0,d,l), 6 = 0.5. 
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